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CO(3— 2) Emission in Spiral Galaxies: Warm Molecular Gas in 

Action? 
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ABSTRACT 

Using the APEX sub- millimeter telescope we have investigated the 12 CO(3— 2) 
emission in five face-on nearby barred spiral galaxies, where three of them are 
high surface brightness galaxies (HSBs) lying at the Freeman limit, and two are 
low surface brightness galaxies (LSBs). We have positive detections for two of 
three HSB spirals and non-detections for the LSBs. For the galaxies with positive 
detection (NGC0521 and PGC070519), the emission is confined to their bulges, 
with velocity dispersions of ~ 90 and ~ 73 km s" 1 and integrated intensities 
of 1.20 and 0.76 K km s _1 , respectively. For the non-detections, the estimated 
upper limit for the integrated intensity is ~ 0.54 K km s _1 . With these figures we 
estimate the H 2 masses as well as the atomic-to-molecular mass ratios. Although 
all the galaxies are barred, we observe 12 CO(3— 2) emission only for galaxies with 
prominent bars. We speculate that bars could dynamically favor the 12 CO(3— 2) 
emission, as a second parameter after surface brightness. Therefore, secular evo- 
lution could play a major role in boosting collisional transitions of molecular gas, 
such as 12 CO(3— 2), especially in LSBs. 

Subject headings: galaxies: ISM; galaxies: spiral; galaxies: stellar content 



1. Introduction 



Low surface brightness galaxies (LSBs) remain among the most intriguing galaxies. De 
fined as galaxies having disk central surface brightness fainter than (B) 22. Qj mag arcsec -2 
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5 Some authors define the limit as 23.0 mag arcsec 2 , for example llmpev &: Bothunl (jl997l ). 
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they are a produc t of low stellar density. LS Bs typically have (1) large amounts of atomic gas 
in the form of HI (jvan der Hulst et al. 1993 ) and (2) low star format ion rates (SFRs). In gen- 
eral t hey have sub-solar metallicity, in agreement with the low SFRs (Ide Blok fc van der Hulst 
19981 ) and consequently the weak production of metals. Also, LSBs have usually large mass 
to light ratios, indicating that disk dynamics is dominated by significant amounts of dark 
matter halos. This is consistent with their flat rotation curves, which in many cases extend 
to several times the optical radii. 

A key questions about LSBs is what physical conditions prevent the gas from forming 
stars. Is this due solely to the fact that the gas is not dense enough to trigger gravitational 
collapse and form stars? Is the amount of molecular gas too low to ensure significant stellar 
formation episodes? How different is the CO-to-H 2 conversion factor X = [N(H 2 )/W(CO)] in 
LSBs compared to the value derived for high surface brightness galaxies (HSBs) in preventing 
their H 2 from being dis covered usin g CO as a tracer? Several studies indicate that A is a 
function of metallicity (Ilsraellll997l ). and thus it should not be a surprise that in LSBs the 
conversion factor reaches higher values than those obtained for HSBs. 

The only way to estimate the amount of molecular gas (H 2 ) in galaxies is to trace the CO 
emission. Only a small number of detections of 12 CO(T— 0) and 12 CO(2— 1) in LSBs have been 
achieved (Oneil. Hofner fc Schinnerer 2000 ; Matthews &: GaolboOl ; O'Neil. Schinnerer fc Hofner 
20031 ; iMatthews et al.ll2005l ). However, a different tracer such as 12 CO(3— 2) e mitted usually 



by warm CO, is necessary to be tter constrain the gas temperature and density (IDumke et al. 
200ll ; iMuraoka fc Kohnoll2006r ). He re we report on the detec tion of 12 CO(3-2) in two HSB 



spirals at the Freeman limiiQ from Galaz et al. ( 2002 . 2006 ). and the nondetection of the 
same transition in another HSB and two LSBs. 

The weak 12 CO(l— 0) emission in LSBs galaxies (see references above) suggests that we 
search for other 12 CO lines better suited for warmer environments. The 12 CO(3— 2) transition 



may b e more easily detected, since it is excited in warm gas (E/k=33.2 K, iMeier et al. 



( I200ll )). which is probably present in LSBs due to the lower metallicity and lack of dust which 
normally shield the UV radiation preventing the excitation of lower CO transitions (10-30 
K). The caveat is whether the UV radiation destroys completely the CO or allows higher 
energy transitions such as 12 CO(3— 2). Because of its higher characteristic temperature and 
critical density (n cr ~ 2 x 10 3 cm -3 ), the 12 CO(3— 2) transition may be more sensitive to 
warm and/or dense gas involved in stellar formation, a key insight that could explain why 
most LSBs have small M# 2 /M#j ratios. 



'Defined as (J.o(B) = 21.65 ± 0.3 mag arcsec 1 ( Freeman 1970h 
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2. The sample 



We have observed five spirals from the sample of iGalaz et al.l (120021 . 120061 ) . All of them 
are face-on spirals, and therefore the APEX beam (of about 18//) samples their bulges. Two 
of the five galaxies are bonna fide LSBs, with disk surface brightnesses fainter than 22.0 mag 
arcsec -2 (UGC02921 and UGC02081), and three are HSB spirals (NGC0521, PGC070519, 
and NGC7589: see below for details). Four of the selected galaxies have measurable near-IR 
emission, and therefore they have a significant population of low-mass and/or evolved stellar 
populations. The criteria to choose galaxies were the following: 



1. face-on orientation - This minimizes the extinction for optical observations, allowing 
us to identify directly the sub-mm line width with the pure gas velocity dispersion. It 
also minimizes the inclination bias on the estimated disk surface brightness. 

2. The HI mass - With the aim of accounting for the atomic gas content as an additional 
variable, we have selected galaxies with different HI masses. The HI masses vary 
between 7 x 10 8 and 51 x 10 8 M m . 



Some information worth mentioning about the selected galaxies follows. 



UGC02081.- Included in the llmpey et al.l ( 119961 ) catalogue of LSBs an d studied by 



Galaz et al.l (12002 . |2006j) (LSB100), it is also in the HIPASS Catalogue dMever et al. 
20041 ). It is at czhei = 2616 km s -1 and has a diameter of 23 kpc. It also exhibits a tiny 



bar in the central region. It has been studied in many optical and near-IR bands, but 
was n ot detected by IR AS. It is not a very massive galaxy in terms of its atomic gas 
mass (IGalaz et al.l 120021 ) . Its disk central surface brightness (B) of 22.4 mag arcsec -2 
locates this galaxy at the LSB re gime. It has a sm all bulge of about 190 pc scale 
length, with colors B — R = 1.11 (IGalaz et al.ll2006l ). The near-IR color of the bulge 
is the same as for the total galaxy, meaning that the bulge stellar populations share 
their properties with those of the disk. 



NGC7589.- Classified as an Sa bv llmpev et al.l (Il996h [LSB473 in Galaz et al. (2002, 
2006)], it has Ho{B) = 21.51 mag arcsec -2 and is therefore an HSB galaxy. Located at 
czheiio = 8938 km s -1 , it has a diameter of 36.12 kpc and has traces of a disk bar. The 
difference between its bulge color (B — R = 2.0) and the disk color (B — R = 1.47) 
shows that this ga laxy has a metal-rich bulge compared to the bulge of other spirals 
( IGalaz et al.ll2006l ). Its near-IR bulge colo r (J — = 0.81) suggests that the bulge is 



more metallic than other bulges in spirals (IGalaz et al.ll2002l ) and larger in size, with a 
scale length of 900 pc (B band). Its HI mass of 51.29 x 1O 8 M makes it ~ 7 times more 
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massive than UGC02081. In spite of its absolute magnitude (Mb = —19.87 mag), the 
galaxy is not detected by IRAS. 



PGC070519.- LSB 463 in iGalaz et al.l (12002 . 120061 ). This SBc galaxy (cz = 5244 km 
s" 1 ) is about 19 kpc diameter. It has a disk central SB /j,q(B) = 21.7 mag arcsec -2 and 
therefore is also an HSB galaxy. This galaxy has an absolute magnitude Mg = —18.63 
and a color B — R = 1.0 in the bulge and B — R = 0.84 in the o yerall galaxy. Thi s 
implies a bulge with larger metallicity compared to that of the disk ( iGalaz et al.ll2002l ). 
for a galaxy with also a large HI mass (37.15 x 1O 8 M ) considering its size. It is not 
detected by IRAS. We note that it also presents a noticeable bar which appears to be 
clearly "melted" with the bulge. 

NGC0521- This is the largest and brightest spiral in our sample (Mb = —20.11). 
Classified as an SBsc(r) (cz = 5018 km s -1 ), has a diameter of ~ 65 kpc, twice that 
of the Milky Way. Its central disk surface brightness is ^q(B) = 21.7 mag arcsec -2 , 
so it is an HSB galaxy. It has visible spiral arms and a prominent bulge with a scale 
length of 720 pc in the R band. The bulge is quite red (B — R = 1.64), with almost 
no difference in color with the overall galaxy (LSB059 in Galaz et al. 2006). The 
near-IR color J — K s = 0.74 suggests a metallic bulge. It has a large amount of HI 
(43.7 x 1O 8 M ). However, as discussed below, this is not a large amount considering 
the remarkable size of the galaxy. It is detected by IRAS in 60 and 100 /im, with fluxes 
of about 0.65 and 3.16 Jy, respectively. It has a noticeable nuclear bar which does not 
however extend too far through the disk. 



UGC02921.- An SAB(s)dm galaxy in llmpev et al.l dl996h (p (B) = 23.6 mag arcsec" 2 ), 
located at czheiio = 3544 km s" 1 ), it is an LSB galaxy. It has a diameter of 21 kpc and 
a large HI mass (21.88 x 10 8 M Q ). It present a tiny bar from which two spiral arms are 
developed. 



3. Observations with APEX 



For APEX observations we use the heterodyne receiver APEX-2A (345 GHz), tunable 
in the frequency range 279 — 381 GHz. The receiver noise temperature (about 60-70 K) 
is fairly constant over the entire tuning range. The telescope beam size at 345 GHz is 
~ 18 arcsec. We observed the five gal axies between 2006 Ju ly and 2007 January, using 
the chopper wheel calibration technique (IKutner fc Ulichl Il98lh . which provides main-beam 
brightness temperature Tmb, after dividing by the main beam efficiency tjmb = 0.73. The 
typical noise system temperature during the observation was T sys ~ 150 — 180 K. The total on 



- 5 - 



source integration time was 2 hr on average, with a velocity resolution per channel of 0.11 km 
s _1 . We tuned the receiver to the corresponding redshifted 12 CO(3— 2) line. Data reduction 
was done with the GILDAS-CLASS package ( |http://www.iram.fr/IRAMFR/GILDAS| ; for 
each galaxy all spectra were added and a linear baseline subtracted. The final spectrum for 
each galaxy was smoothed to obtain a velocity resolution of ~ 16 km s -1 and an rms noise 
temperature T rms = 5 mK. Therefore, the typical noise temperature per channel is about 
5 x ^16/0.11 ~ 60 mK. 

Figure 1 shows final smoothed and summed spectra. Only for galaxies NGC0521 and 
PGC070519 we detected 12 CO(3— 2) emission. A Gaussian fit was applied to each spectrum 
and the parameters of the fit are summarized in Table [TJ Detections are defined for signals 
above 3a, where a ~ 5 mK is the rms noise temperature. Galaxies NGC7589, UGC02081, 
and UGC02921 do not present 12 CO(3— 2) emissions larger than the detection threshold. 



4. Analysis and discussion 



Driven by the uncertainty in t he CO-to-H? conversion facto r for s pirals in general, we 
use an approach similar to that of lO'Neil. Bothun fc Schombertl (120001 ) to compute the H 2 
masses for NGC0521 and PGC070519 and the corresponding upper limits for NGC7589, 
UGC02081, and UGC02921. Therefore, we assume a 12 CO(1-0) to 12 CO(3-2) ratio of 1. 
This is a reaso nable approa c h con sidering, for example, the ratio of about 1.2 ob t ained , 



on average, by iMeier et al.l (120011 ) for a sample of dwarf galaxies. iDumke et al.l (120011 ) 
obtained similar figures, with a (3-2)/(l-0) ~ 1.3 for the centers of spirals. CO velocity 
dispersions are ~ 89 km s _1 for NGC0521 and 73 km s _1 for PGC070519. Since these 
galaxies are face-on, the velocity width would correspond to the intrinsic velocity dispersion 
of the 12 CO(3— 2) emission. In both cases, and given the similar sizes of the corresponding 
bulges and the antenna beam, the velocity dispersions would simply correspond to that 
of the CO embedded in the respective bulges. Note that the velocity dispersions are ~ 
5 — 10 times larger than those typically observed for HI in the central regions of spirals 
(Ivan der Kruit fc Shostaklll982l ). We suspect that this large width is due to velocity fields 
induced by the bars observed in NGC05 21 and PGC070519 (see Figure 1). Such behavior 
was also observed by lDumke et al.l (120011 ). who obtained similar values for AV^. In Table[T]we 
indicate the 12 CO(3— 2) line velocity width and also the HI velocity dispersion. We note the 
high value also for W(HI) in NGC0521, strongly indicating the presence of a strong velocity 
field. Note that all the other galaxies present large velocity dispersions in HI, suggesting 
therefore a likely large value for the velocity dispersion of their molecular gas. 



From the Gaussian fits we determine the integrated intensity \qo — J T MB dv, and 
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Fig. 1. — Images and spectra of galaxies observed with APEX. Images are are from the SDSS 
(B filter). Spectra are centered in the 12 CO(3— 2) rest emission. The velocity resolution was 
smoothed to about 16 km s _1 for each spectrum. The solid line for NGC0521 and PGC070519 
spectra indicate the corresponding Gaussian fits to the line emission. The circle in each image 
represents the mean antenna beam size of about 18 arcsec. 
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the intensity at the pe ak of the emission (in mK, see Table 



Bregman &: Hoga (119881 ) and the formula of ISanders et al.l (119861 ) 



TJ). Using the method of 



M T {H 2 ) = 5.82[vr/4]^j co , (1) 

where <4(pc) is the telescope beam diameter at the distance of the source and Ico is the total 
CO line integrated intensity (K km s _1 ), we estimate the H 2 mass for galaxies with positive 
12 CO(3-2) detections (NGC0521 and PGC070519). For galaxies with negative detections 
(UGC02081, UGC02921 and NGC7589) we estimate upper limits using 

Ico < 3T MB Av H i/ \fn, (2) 

in K km s -1 . We assume that AVhi ~< AVi2co(j=3-2) >, the average velocity dispersion 
between NGC0521 and PGC070519, that is, 80 km s -1 . We think this value is more realistic 
than just using an arbitrarily smaller value, since UGC02081, UGC02921, and NGC7589 
also present bars. Thus y/n = a/80/16 = \/5, the number of channels used in the smoothed 
spectra. Therefore 

3 

Ico < —7=T MB Av co . (3) 
V5 

For non-detections (UGC02921, UGC02081 and NGC7589), we use T MB = a rms ~ 5 mK for 
all the corresponding spectra. Thus, Ico < 0.54 K km s _1 . Using the corresponding beam 
diameter in pc for all galaxies, as seen at the distance of each galaxy, we obtain estimated 
H 2 masses for NGC0521 and PGC070519, and upper limits for UGC02921, UGC02081, and 
NGC7589 (see eq. [1] and Tabled). It is worth noting that we do detect 12 CO(3-2) for 
NGC0521 and PGC070519, but do not detect it for NGC7589, which is 0.2 mag arcsec" 2 
brighter than these two galaxies, suggesting that factors other than surface brightness must 
play a key role in the 12 CO(3— 2) emission. 

Tabled] presents H2 masses. Note that these values are computed assuming that molec- 
ular gas is uniformly distributed through galaxies, which we know is not the case, and in all 
galaxies, the beam was pointed to the bulge. As shown in Table [TJ H 2 masses derived from 
detections are ~ 10 8 M . The estimated upper limit for the molecular mass in the beam 
size of galaxies with no detections is ~ 3 - 28 x 10 7 M . With these values we are able to 
compute the molecular-to-atomic gas mass fractions for each galaxy (again, only upper limit 
estimates for galaxies with no detections). The results agree with the picture that LSBs lack 
significant amounts of molecular gas and, in general, with small molecular-to-atomic mass 
ratios (below 0.08). For one case (NGC0521) the molecular-to-atomic gas fraction appear 
as large as 0.9. Although surprising, this result could be anticipated given its small HI 
mass considering its large size. In other terms, the molecular gas content for NGC0521 is 
comparable to its atomic one. Note that we do not detect 12 CO(3— 2) for three but instead 
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for two HSBs and do not detect any CO for both of the LSBs. How could one explain this 
puzzling picture? Although the central disk SB of these galaxies seems to play a role in the 
CO emission, we suspect that other factors are key in allowing such an emission. Looking in 
detail at the structure of the galaxies, we note that NGC0521 and PGC070519 have much 
more prominent bars really competing with the bulge emission, compared to the other three 
galaxies (see Fig. 1 and also figures in Galaz et al. 2006), including NGC7589, which is 
0.3 mag arcsec" 2 brighter than NGC0521 and PGC070519. We suggest that strong velocity 
fields could be responsible for the molecular emission in LSBs. 



5. Conclusions 

We detect 12 CO(3— 2) in 2 of 5 nearby spirals, using the APEX submillimeter telescope. 
The emission is detected for NGC0521 and PGC070519, both HSB galaxies lying at the 
Freeman limit. The other galaxies, with no detections, share similar morphology and ori- 
entation, but two are LSBs (UGC02921 and UGC02081) and one an HSB (NGC7589). All 
of them are face-on, making the internal extinction negligible. The measured main-beam 
temperature CO fluxes are 1.20 K km s" 1 for NGC0521 and 0.76 K km s" 1 for PGC070519. 
The remaining galaxies have fluxes below 3a, and thus we are able only to estimate upper 
limits for their 12 CO(3-2) fluxes (< 0.54 K km s^ 1 ). 

Measured velocity dispersions for NGC0521 and PGC070519 are 89 and 73 km s _1 , 
respectively, ~ 5 — 10 times typical values obtained by other authors for bulges in spirals 



(van der Kruit & Shostak 


1982; 


Dumke et al. 


2001) 



velocity field is dominated by the bar kinematics (IDumke et al.ll200ll ). Bars are observed 
actually in all five galaxies, but those of NGC0521 and PGC070519 appear as the more 
prominent ones. 

We compute the total H 2 mass in the main beam, obtaining 2 x 10 8 M Q and 1.4 x 1O 8 M 
for NGC0521 and PGC070519, respectively. The corresponding molecular-to-atomic gas 
fraction is about 0.9 and 0.08. The high value for NGC0521 is probably due to its low HI 
density, and to the assumption that the HI is uniformly distributed over the whole galaxy 
optical size. In fact, in spirals one expects that the HI is mostly distributed along the disk. 

Overall, we have shown that 12 CO(3— 2) emission could be intense for some galaxies at 
the SB "Freeman limit." Moreover, such an emission could be "dynamically boosted" by 
bars, helping to warm the CO and allowing the 12 CO(3— 2) transition at 33 K above the 
ground level. Although we have no detections of 12 CO(3— 2) for the LSBs in this sample, 
we speculate that many LSBs could present 12 CO(3— 2) emission given some dynamical 
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conditions that warm the gas (or make it denser). We speculate that many LSBs with no 
12 CO(l— 0) or 12 CO(2— 1) emission could have instead 12 CO(3-2) emission thanks to the 
poor UV shielding favored by the low metallicity and low dust content, allowing a higher 
gas temperature. This emission, which in principle could be small, may be largely amplified 
by secular processes which warm the gas. 

G.G. acknowledges support from FONDECYT 1040359. P.C. was funded by the CONICYT- 
ALMA Fund, project 31050003. L.B. and M.R. acknowledge support from the Center for 
Astrophysics FONDAP 15010003. We are grateful to the anonymous referee for helping to 
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Table 1: Gaussian fit parameters for 12 CO(3— 2) emission in NGC0521 and PGC070519, and upper limits for non- 
detections, 



Name 


Tmb 




J T MB dv 




M# 2 (18 arcsec) 


M Hl (18 arcsec) 


M H JM HI 


Width 


W(HI) 




(mK) 


(mK) 


(K km s" 1 ) 


(xlO 8 M 


) (xlO 8 M ) 


(xlO 8 M ) 




(km s~ r 


) (km s _1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


10) 


NGC0521 


12.7 


2.9 


1.20 ±0.10 


43.7 


2.01 


2.2 


0.91 


88.7 


244 


PGC070519 


9.8 


2.5 


0.76 ±0.09 


37.1 


1.39 


16.5 


0.08 


72.8 


64 


NGC7589 


< 5 


-5 


< 0.54 


51.3 


< 2.77 


26.6 


< 0.10 


~ 80 


183 U 


UGC02081 


< 5 


~ 5 


< 0.54 


7.76 


< 0.25 


1.1 


< 0.23 


~ 80 


179 i 


UGC02921 


< 5 


~ 5 


< 0.54 


21.8 


< 0.45 


13.4 


< 0.03 


~ 80 


168 


Notes: (1) From |Gali 


iz et al. (2002, 20061) and 


the NED. (2) Calit 


>ratcd peak main- 


beam brightness temperature. 


(3) The rms estimates are 


from the smoothed 


data with 16 


km s 1 



resolution. (4) Integrated intensity. (5) HI mass fromjGalaz ct al. (200^) and the NED. (6) Estimated H2 mass enclosed by the APEX main-beam of 18 arcsec. Values for last 3 galaxies 



arc upper limits. (7) HI mass interpolated to the APEX main-beam. (8) H2 to HI mass ratio. (9) Velocity width of the line. For the last 3 galaxies the velocity width is estimated as the 
average value obtained for NGC0521 and PGC070519. (10) Velocity width for the HI line, obtained from HIPASS jMeyer et"al~||2004| ) and HYPERLEDA jPaturel et al.||2003l ), clipped at 
20% peak flux density. For PGC070519, and NGC7589 the indicated value correspond to the mean homogenized maximum rotation velocity uncorrected for inclination. 



